Dynamics of entanglement in quantum computers with imperfections 
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The dynamics of the pairwise entanglement in a qubit lattice in the presence of static imperfections 
exhibits different regimes. We show that there is a transition from a perturbative region, where the 
entanglement is stable against imperfections, to the ergodic regime, in which a pair of qubits becomes 
entangled with the rest of the lattice and the pairwise entanglement drops to zero. The transition 
is almost independent of the size of the quantum computer. We consider both the case of an initial 
maximally entangled and separable state. In this last case there is a broad crossover region in which 
the computer imperfections can be used to create a significant amount of pairwise entanglement. 
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Any practical implementation of a quantum computer 
will have to face errors, due to the inevitable coupling 
to the environment P, Q or to device imperfections . 
The effect of the environment is to introduce decoher- 
ence which sets the time scale over which quantum com- 
putation is no longer possible. The presence of static 
imperfections, although not leading to any decoherence, 
may be also decremcntal for the implementation of any 
quantum computational task. A small inaccuracy in the 
coupling constants induces errors in gating or an un- 
wanted time evolution in the case in which the Hamil- 
tonian cannot be switched exactly to zero. If the im- 
perfection strength increases, new phenomena occur and 
above a certain threshold the core of the computer can 
even "melt" due to the setting in of chaotic behavior Q . 
Previous investigations of the stability of quantum infor- 
mation processing in the presence of such effects mainly 
studied the fidelity of the quantum evolution as an in- 
dicator of the quality of the computation 0, 0, 0|- In 
particular, in Ref. l3| the fidelity was used to measure 
the stability of the quantum memory, that is of a state 
loaded on a quantum computer with imperfections. A 
more complete characterization of the stability of a quan- 
tum computation requires a deeper investigation of the 
state of the system. Fidelity is one of the ways to char- 
acterize it. In this Letter we discuss the behavior of en- 
tanglement on approaching the transition to quantum 
chaos. Entanglement is not only one of the most intrigu- 
ing features predicted by the quantum theory but also 
a fundamental resource for quantum computation and 
communication 0. Therefore studies of the stability of 
entanglement under decoherence and imperfection effects 
are, in our opinion, of interest. 

We model the quantum computer as a lattice of in- 
teracting spins (qubits) where, due to the unavoidable 
presence of imperfections, the spacing between the up 
and down states (external field) and the couplings be- 
tween the qubits (exchange interactions) are both ran- 
dom. The entanglement properties in spin systems has 



recently attracted attention and we refer to 7] for a more 
detailed introduction. In this work we consider, for the 
first time, the effect of disorder in the couplings. 

We consider n qubits on a two-dimensional lattice, de- 
scribed by the Hamiltonian 

(1) 

where the erf {a = x, y, z) are the Pauli matrices for the 
qubit i and the second sum runs over nearest-neighbor 
qubit pairs on the lattice. The energy spacing between 
the up and down states of a qubit is given by = 
Ao -f 5i, where the 6iS are uniformly distributed in the 
interval [—(5/2,(5/2]. The parameter 5 gives the width of 
the Ai (single qubit energy) distribution around the av- 
erage value Ao . The couplings Jy are also uniformly dis- 
tributed in the interval [—J, J]. We consider < 7 < 1. 
For 7 = 1 Eq. ^ reduces to the disordered Ising model, 
while for 7 = it gives the disordered XY model. For 
7=1 the model defined by Eq. has been proposed by 
Georgeot and Shepelyansky ^] to describe the hardware 
of a quantum computer, in which system imperfections 
generate unwanted interqubit couplings and energy fluc- 
tuations. For J, (5 = the spectrum of the Hamiltonian 
is composed of n + 1 degenerate levels, and the interlevel 
spacing is 2Ao, which corresponds to the energy required 
to flip a single qubit. We study the case (5, J ^ Ao, m 
which the degeneracies are resolved and the spectrum is 
composed by n+1 bands. Since d,J^ Aq, the interband 
couplings in can be neglected and every single band 
can be studied separately. We concentrate our studies 
on the central band with zero magnetization (^ a^ = 0). 
We assume free boundary conditions and set the energy 
scale Ao — 1. 

In this Letter, we investigate the evolution in time of 
the entanglement between two nearest neighbor qubits in 
the lattice. We confine our interest to the entanglement 
of formation between two spins in the lattice. The 
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propagation of entanglement in disordered-free systems 
has been analyzed in Ref. 0- In the same spirit we 
study the entanglement evolution for an initially maxi- 
mally entangled pair of qubits: 

l^BiO)) - -^(|01) + |10)) ® 10101 ... 01) . (2) 

We will also analyze the generation of the entanglement, 
due to interqubit coupling, of two spins which are initially 
in a separable state. In this case we start from the wave 
vector 

|*s(0)> = |01) ® 10101... 01) . (3) 

In both cases, the initial state, Eqs. I|2I3|I . evolves ac- 
cording to the Hamiltonian defined in Eq. (Q. Wc 
compute the reduced density matrix pi2{t), obtained 
from the wave function |vE'(i)) after all the spins except 
those at sites 1 and 2 have been traced out: pi2{t) — 
Tr3,...,„|*(t))(^'(t)|. We choose qubits 1 and 2 to be 
nearest neighbors on the border of the lattice (their po- 
sition in the lattice is not crucial for the results here 
reported). Their state at time t, described by the den- 
sity matrix /0i2(i), is mixed since qubits 1 and 2 become 
entangled with the rest of the lattice. We quantify the 
entanglement by means of the concurrence j3| , defined as 
C — max{Ai — A2 — A3 — A4, 0}, where the Ai's are the 
square roots of the eigenvalues of the matrix R = P12P12, 
in descending order; the spin-flipped density matrix is 
defined by pi2 — {al (g) 0'2)Pi2('''i ^ '^2) i^^ this defi- 
nition the standard basis {|00), |01), |10), |11)} must be 
used) . Another dynamical quantity widely used to char- 
acterize the stability of a state under perturbations is 
the fidelity 3]. It is defined as f{t) = |(*(0)|*(i))p, 
where \^{t)) = ex.p{-iHt)\->p{0)) (we set h = 1). For 
weak interqubit coupling f{t) is close to 1 for all times. 
For strong coupling we enter the so-called quantum chaos 
regime, and the fidelity essentially drops to zero Q- In 
the following we contrast the behavior of the fidelity with 
that of the concurrence. In addition to the different 
regimes found in the analysis of the fidelity, the concur- 
rence gives additional information related to the choice 
of the initial state. 

Initially entangled state - We first consider the time 
evolution of C{t) starting from the (maximally entan- 
gled) initial state |^'b(0)). If the coupling strength J 
is weak, then the time evolution of C{t) can be ob- 
tained by means of perturbation theory in the coupling 
J. Only the quantum register states directly coupled 
to the initial state are involved (the quantum register 
states are the eigenstates of Hamiltonian at J = 0: 
|a) = |ai, Q!2, an), with a,; = 0, 1). Since the exchange 
interaction has a two-body nature, the matrix elements 
are zero when the state binary strings differ in more than 
two digits. Furthermore, only nearest neighbor qubits in- 
teract, and therefore to first order in perturbation theory 
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FIG. 1: Concurrence time scale tc as a function of J, for 
7 = 1, (5 = 0.2, initial state \'^b{0)), and n — 4 (triangles up), 
6 (triangles left), 8 (triangles down), 10 (circles), 12 (squares), 
and 14 (diamonds) qubits. Stars represent n = 10 and 7 = 
0. The full and the dashed lines are proportional to 
and J~^, respectively. Inset: Concurrence as a function of 
time for n = 10 and 7 = 1. From top to bottom: J = 
2, 5 X 10~^; 1, 2, 3 X 10~^. To smooth statistical fluctuations, 
here and in the following figures data are averaged over A'^,- 
different random configurations of the imperfections Si and 
Jij : Nr^50 for n = 4, 6, 8, 10 and Nr = 30, 20 for n = 12, 14, 
respectively. 

only order of n quantum register states are involved. It 
is important to note that the coupling to a state which 
differs from the initial state |^b(0)) in the polarization 
state of two qubits other that 1 and 2 does not affect the 
pairwise entanglement between qubits 1 and 2. To first 
order in the perturbation theory, these two qubits become 
entangled only with nearest neighbor qubits. Therefore 
the concurrence can be described in terms of superposi- 
tion between a small number of quantum register states, 
independently of the total number n of qubits. In the 
perturbative regime, the frequencies of these oscillations 
are ^ S, around an average value ~ 1 — (J/ 5)"^. Perturba- 
tion theory breaks down when the typical interaction ma- 
trix element J between directly coupled states becomes 
of the order of their energy separation AE S/n. This 
follows from the fact that two quantum register states 
are directly coupled when two qubits are flipped. This 
changes the energy by ~ (5, and there are the order of n 
states directly coupled inside this energy interval. There- 
fore the perturbative regime is limited to J < Jp ^ S/n 
0. Above this threshold, it is known i.'l^ that the fi- 
delity decays exponentially: /(t) « exp(— F/t), with the 
rate F^^ ~ J^Pf determined by the Fermi golden rule, 
where Pf — 1/AE ~ n/6 is the density of directly cou- 
pled states. The transitions involving qubits other than 1 
and 2 do not change the concurrence. Therefore the den- 
sity of directly coupled states relevant for this quantity 
is given by pc ^ 1/6. Hence we expect the concurrence 
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FIG. 2: Concurrence saturation value Coo as a function of 
J with same parameter values and meaning of the symbols 
as in Fig. 1. The dashed curve represents the fit Coo (J) ~ 
exp(-A(J - Jo)), with ^ « 20 and Jo « 1.0 x 10"^. 




FIG. 3: Concurrence C{t) for 7 = 1 (full curves) and 7 = 
(dashed gray curve), 5 — 0.2, initial state |'I's(0)), n — 10, 
and, from bottom to top, J = 1, 2 x 10"*; 1, 2 x 10~^; 1, 2 x 
10"^; 1, 2 X 10"^ The straight lines give the time dependence 
C{t) — Jt predicted by the perturbation theory. 



to decay exponentially, C{t) ~ exp(— Fct), with the rate 
Fc ~ J^Pc ~ J^/(5. This expectation is borne out by the 
numerical data of Fig. 1, which show that the concurrence 
time scale , defined by the condition C(if ) = 0.96, is 
inversely proportional to (the result does not depend 
on the value chosen). In the regime of strong coupling, 
J > 5, the time evolution of the fidelity is given by |3l.ll0l| 
f{t) « Y[i<j cos{Jijt), where the product runs over near- 
est neighbor qubits. Thus the fidelity decay is Gaussian: 
/(<) « exp(- Y.i<j Jfji^) ~ exp(- J^ni^). In computing 
the concurrence decay, we must again take into account 
that the rotations of qubits different from 1 and 2 do not 
change their pairwise entanglement. Therefore we ex- 
pect C{t) ^ exp(— J^t^). This gives a concurrence time 
scale oc J~^, in good agreement with the numerical 
data of Fig. 1. Since in the Fermi golden rule regime 
the concurrence time scale ^ ^/J'^, while in the er- 
godic regime ^ 1/J, the transition between these two 
regimes takes place at J = ~ 5 (we checked numeri- 
cally that Je oc S). The data of Fig. 1, taken at different 
number of qubits, seem to collapse on a single curve, in 
agreement with our theoretical expectations. 

It is also interesting to study the saturation value of 
the concurrence. Coo = (limt^oo C(i')dt')/r. In or- 
der to smoothen the fluctuations, we average our numer- 
ical data over both different random configurations of 
di and Jij and a large enough time interval r (we take 
T equal to 1/10 of the total integration time t, with t 
sufficiently long to observe the saturation of C{t)). In 
Fig. 2, we study Coo as a function of the coupHng strength 
J, starting from the entangled state |^'b(0)). Again we 
can distinguish three regimes. In the perturbative regime 
J < Jp ^ 6/n the saturation value is close to one (more 
precisely, 1 — Coo oc (J/S)"^). In the ergodic regime at 
J > Je S the concurrence Coo is equal to zero. Quan- 



tum chaos is characterized by the ergodicity of the eigen- 
functions {i = 1,...,N = 2") of Hamiltonian JQ. 
This means that, by expanding over the quantum 
register states, = X^aCala), the coefficient Cq are 
randomly fluctuating and have amplitudes |cq.| ^ 1/^/N. 
In this regime, the complexity of the eigenstates, char- 
acterized by the entropy Si = — X]q I'^aP log |cap, is 
maximal, that is Si ~ n. We note that indeed the 
maximum entropy criterion iSj « n gives the ergodic- 
ity threshold Je ^ S 0|. Due to quantum ergodicity, 
for times t >> if ~ 1/J the wave function |\l/(i)) is 
a random superposition of the quantum register states: 
= J2a^a{t)\a), where the coefficients aa{t) have 
amplitudes ^ and random phases. The diagonal 

elements of the reduced density matrix pi2 are given by 

[Pl2(t)] ai .OL2 ',oli -,(^2 = Ea3,...,a„ |aai....,a„(t)|'. The valuc 
of these elements is « 1 /4, since they are the sum of iV/4 
positive terms, and the value of each term is ~ 1/A''. 
The off diagonal matrix elements of pi2 are instead given 
by the sum of N/A terms with random signs, and there- 
fore their value is ^ For such a reduced den- 
sity matrix the concurrence is equal to zero. Physically 
this means that in the ergodic regime the entanglement 
is multipartite and shared between all the qubits, and 
therefore the pairwise entanglement drops to zero. In 
this regime the coupling to the rest of the system acts as 
a dephasing channel, and the reduced density matrix be- 
comes essentially diagonal, thus destroying the pairwise 
entanglement. 

Initially separable state - Since the Hamiltonian in 
Eq. can create entanglement due to the exchange 
term between neighboring qubits 0, it is interesting to 
consider also separable states as initial states. In the 
perturbative regime J < Jp, the concurrence oscillates 
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FIG. 4: Concurrence saturation value Coo as a function of J 
for 7 = (black symbols) and 7 = 1 (gray stars), 5 = 0.2, 
initial state \^ s) (filled symbols). Symbols have the same 
meaning as in Fig. 1. The straight line is proportional to J. 
The dashed line is the same as in Fig. 2. 

with frequency | (Si — 1^2 1 between and a maximum value 
2|Ji2|/|<5i - 52\. At times t < l/\5i - 52\, the pertur- 
bation theory gives C[t) « 2|Ji2|t, in agreement with 
the numerical data of Fig. 3 (after disorder averaging we 
get C(t) = Jt). In the ergodic regime J > Je, a max- 
imum concurrence of the order of one is generated after 
a time ~ 1/J. After this time the entanglement of the 
pair rapidly drops to zero. In the intermediate regime 
the concurrence reaches its maximum and then saturates 
to a value which is nearly independent of J. Therefore 
the generation of pairwise entanglement is an interesting 
characteristic of this region. 

The saturation value Coo as a function of the cou- 
pling J is shown in Fig. 4. In the perturbative regime 
Coo oc J/5 whereas in the ergodic regime Coo =0. It 
is interesting to notice that there is a broad region in 
which the interactions create a significant amount of en- 
tanglement, quantified by a concurrence Coo ~ 0.2 — 0.3. 
We also note that in the ergodic regime J > Je the sat- 
uration value Coo = does not depend on the initial 
conditions, as expected. On the contrary, when J < Je 
the concurrence saturation value depends on the initial 
state vector, implying the absence of ergodicity. 

As it can be seen from our data shown in Figs. 1-4, we 
did not found any significant change in the dynamics of 
entanglement as a function of 7. Even though a further 
symmetry is added to the system when 7 = (the in- 
variance with respect to rotations about the z-axis) , this 
has little influence on entanglement. 

In summary, we have studied the dynamics of the 
pairwise entanglement in a qubit lattice in the presence 
of static imperfections, and characterized three different 
regimes: (i) the perturbative regime, in which the entan- 
glement is stable against imperfections, (ii) the crossover 



regime, in which the imperfections degrade the concur- 
rence of an initially entangled pair but can also drive a 
significant entanglement generation, and (iii) the ergodic 
regime, in which a pair of qubits becomes entangled with 
the rest of the lattice and therefore the concurrence of the 
pair drops to zero. We stress two important points of our 
findings from the point of view of quantum computation. 
First of all, the pairwise entanglement is destroyed above 
a coupling strength which is independent of the size of the 
quantum computer. Moreover, there is a broad crossover 
region in which the computer imperfections can be used 
to create a significant amount of pairwise entanglement. 
Finally we want to emphasize that our analysis leads to 
the conclusion that the spin system behaves as a dephas- 
ing environment only in the ergodic regime. In the inte- 
grable limit, even an infinite system is unable to dephase 
completely the initial state. We believe that this obser- 
vation may be relevant in understanding the behavior of 
spin or more complex baths. 
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